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A truncated determinant algorithm is used to study the physical effects of the quark eigenmodes associated 
with eigenvalues below 420 MeV. This initial high statistics study focuses on coarse (6 4 ) lattices (with 0(a 2 ) 
improved gauge action), light internal quark masses and large physical volumes. Three features of full QCD are 
examined: topological charge distributions, string breaking as observed in the static energy and the eta prime 
mass. 



1. TDA Simulations on Coarse Lattices 

In the truncated determinant approach (TDA) 
to full QCD, the quark determinant, T>(A) = 
Det(H) = Det (75(^4.) - m)) is split-up gauge 
invariantly into an infrared part and an ultravio- 
let part (I . 



V{A)=V IR {A)V UV {A) 



(1) 



The ultraviolet part, 'DjjVi can be accurately fit 
by a linear combination of a small number of Wil- 
son loops 0]. The infrared part T>ir(A) is defined 
as the product of the lowest N\ positive and neg- 
ative eigenvalues of H, with |A$| < Aqcd- The 
eigenvalues A a of H are gauge invariant and mea- 
sure quark off-shellness. The cutoff (for the sep- 
aration in Eq. [l]) is tuned to include as much 
as possible of the important low-energy chiral 
physics of the unquenched theory while leaving 
the fluctuations of In T>j r of order unity after each 
sweep updating all links with the pure gauge ac- 
tion (assuring a tolerable acceptance rate for the 
accept/reject stage) [jl|,{§. This procedure works 
well even for kappa values arbitrarily close to 
kappa critical. 

Initial studies using TDA focus on the qualita- 
tive physical effects of the inclusion of the infrared 
quark eigenmodes. For this purpose, coarse lat- 
tices with large physical volumes are appropriate. 
Here we report the study of 6 4 lattices with an 
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0(a 2 ) improved gauge action 
S G = /?pi aq 51 ^ReTr(l - U 

plaq 

+/3trt^7:ReTr(l-Ut r t) 



plaq J 



(2) 



with /3 p i aq = 3.7 and A rt //3 p iaq = 1-04. This 
0(a 2 ) improved gauge action was adjusted (as in 
Ref. Q]) to have the lattice scale a ss 0.4 fm. 

We used naive Wilson fermions with light quark 
masses (k = 0.2050, 0.2044 and 0.2038); and keep 
the lowest 840 eigenvalues of V, i.e. A cut « 420 
MeV (1/a = .5 GeV). No terms modelling V uv 
were included in the gauge action. 

With k = .2050, approximately 15,000 full 
steps were required for the lattice configurations 
to equilibrate (reflecting critical slowing down - 
a few hundred suffice on small physical volumes). 
The equilibrated behaviour of T>m is shown for 
50,000 sweeps in Figure [|. 

In the remaining sections we consider some 
physical quantities that differentiate between 
quenched and full QCD. 



2. Topolopy 

The topological charge, Q, can be expressed 
in terms of the eigenvalues of the Wilson-Dirac 
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Figure 1. Equilibrated quark determinant in 
TDA simulations on a coarse 6 4 lattice, (k = 
.2050) The lowest 840 eigenvalues included in 
In Det H. 



Figure 2. Topological charge distribution for 
k = .2050 on a coarse 6 4 lattice (dots). Expected 
chiral behaviour (diamonds) is also shown. 



operator. 

1 N 1 

i— l 

This sum is quickly saturated by the low eigen- 
values. 

In full QCD configurations very small eigenval- 
ues of H are suppressed by the quark determi- 
nant factor. In particular, non-zero topological 
charges must be suppressed in the chiral limit 
(m q — > 0). Furthermore, the functional depen- 
dence of the topological charge distribution, Pq, 
on the light quark mass m q is predicted by the 
chiral analysis of Leutwyler and Smilga^j. 

P Q = I Q {xf - I Q+1 (x)I Q -i{x) (4) 

where x = l/2Vf^m^ — Vm q < iptp >, Iq are 
modified Bessel Functions of order Q and V is the 
total space-time volume. 

The qualitative agreement with the expected 
behaviour of Eq. |] has been reported previ- 
ously [Q. Using the value of i = 3.8 deter- 
mined from the lattice calculations (m^ = .39 
and f n = .20 for k = .2050), excellent quantita- 
tive agreement between theory and lattice calcu- 
lations is shown in Figure ||. The data from a 
single kappa value (e.g. k = .2050) determines 



k c . Equating the calculated and predicted be- 
haviour we find the factor 1 — k/k c = 0.0082, i.e. 
k c = .2067. 

3. Static Energy 

The static energy V(R) is defined as: 

V(R)= hmln^r- 1 ) (5) 
t^oo W(R,T) 

In Coulomb gauge, W(R, T) is the Wilson line 
correlator and V(R) is the energy of the lowest 
physical state coupling the static quark-antiquark 
state (with separation R) to the vacuum. As T — > 
oo, V(R) is approached from above. 

Light quark masses implies long autocorrela- 
tions in physical quantities. However, with a bin 
size of 250 sweeps, the values of W((R),T) are 
decorrelated for all relevant R and T||. 

The results for the static energy (n — 0.2050) 
are shown in Figure [| Clear evidence for string 
breaking is seen for T > 3. This strengthens pre- 
vious indications for string breaking ||. 

We can also extract the Sommer scale 
(r%[V(r) - V(r - d)]/\d\ = 1.65) from our de- 
termination of the static energy. Very roughly, 
r w 1.25 and 1/a w 500 MeV. 
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Figure 3. Static energy measured using time 
slices: [0:1] (x's), [1:2] (squares), [2:3] (circles) 
and [3:4] (diamonds), (« = 0.2050) The heavy- 
light pair threshold is denoted by the dashed line. 



Figure 4. The total rj propagator for k = .2050 
on the 6 4 lattices. The valence term (circles), 
hairpins term (squares) and total propagator (di- 
amonds) are shown separately. 



4. Spectrum of Light Hadrons 

Masses for the light hadrons can be determined 
for the various kappa values. For n = .2050 
(highest statistics) we obtain: m T = 0.391(6), 
m p = 1.008(8) and m P = 1.935(89). Significant 
scale violation is seen with naive Wilson fermions. 
For example, the scale 1/a determined from the 
proton mass is 485 MeV while the rho mass gives 
764 MeV. 

The relation between the axial U(l) anomaly 
and the rj mass is well understood in full QCD. 
For two light quarks (Nf — 2), mi = m 2 + m\ 
where = 2Nfx/ f% and the topological suscep- 
tibility is V\ =< Q 2 > quenched. The full t] prop- 
agator is the sum of a connected (valence quark) 
term and a disconnected (hairpin) term. Thus, in 
the continuum, the momentum space full propa- 
gator can be written: 



(p 



rn 

,2 



)- l = 



2 



2 +ml)- 1 (6) 

These separate terms and their sum are shown 
in Figure |J for k = 0.2050. The cancellation be- 
tween the valence and hairpin terms in the full 
propagator is evident. We obtain m v ' = 1.16(12). 
As expected the rj mass is much heavier than the 
pion mass. 



5. Summary and Outlook 

Three results are clear from this study: (1) 
The topological charge distribution behaves as 
expected in full QCD. (2) Static energy shows 
clear evidence for string breaking for T > 1 fm. 
(3) Light hadron spectrum results on a 6 4 coarse 
lattice look reasonable. In particular, m, ; <a = 
1.16(12) > m^a = .391(6). 

The full parametric studies needed to optimize 
with TDA method will require many runs on 
small sets of CPU's. This might be done effec- 
tively by a seti@home approach |J. 
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